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We present our first nucleosynthesis results from a numerical simulation of the thermonuclear 
disruption of a static cold Chandrasekhar-mass 12 C/ 16 white dwarf. The two-dimensional sim- 
ulation was performed with an adaptive-mesh Eulerian hydrodynamics code, Flash, that uses as 
a flame capturing scheme the evolution of a passive scaler. To compute the isotopic yields and 
their velocity distribution, 10,000 massless tracer particles are embedded in the star. The parti- 
cles are advected along streamlines and provide a Lagrangian description of the explosion. We 
briefly describe our verification tests and preliminary results from post-processing the particle 
trajectories in (p, T) with a modest (214 isotopes) reaction network. 

1. Introduction 

Type la supernovae are bright stellar explosions that are distinguished from other supernovae 
by a lack of hydrogen in their spectra. A widely accepted progenitor for these events is a 
close binary system in which the primary component is a degenerate C-0 rich white dwarf 
that accretes from a low-mass companion. The explosion is believed to ignite in the center of 
the white dwarf when its center becomes sufficiently dense and hot. For a review of different 
explosion models, see Q and references therein. 

Constraints on the nature of the explosion come from studies of the supernova light curve and 
the composition of the ejecta. Multi-dimensional hydrodynamical simulations cannot include 
more than a few isotopes and are typically Eulerian. Here we report on our implementation 
of tracer particles into a simulation of a deflagrating C-0 white dwarf. These particles are 
advected along streamlines and provide a Lagrangian description of the explosion that can be 
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post-processed with a realistic reaction network. We discuss our verification tests, a 2-D defla- 
gration, and the resulting nucleo synthetic yields. 

2. Simulation with Tracer Particles 

The simulations of the deflagrating white dwarf were performed with the Flash code 0. 
Flash solves the Euler equations for compressible flow and the Poisson equation for self- 
gravity on an adaptive mesh. These simulations employ a flame capturing scheme (see [3 1) that 
advances the flame by evolving a passive scalar variable with an advection-reaction-diffusion 
equation. The flame model assumes that the flame propagates at a given speed v = max(vi am , v tur b), 
where v\ am is the laminar flame speed [4], and v tur b is a turbulent flame speed based on the as- 
sumption that the turbulent burning on macroscopic scales is driven by Rayleigh-Taylor insta- 
bility [3 1 . The scalar variable is used to model the burning from 12 C to 24 Mg. Further burning to 
Si-group and relaxation to NSE (here taken to be 56 Ni) are modeled as two consecutive stages 
of burning following the burn from C to Mg. 

In addition to the flame model, the simulation of the deflagration had 10,000 tracer particles 
added to track the evolution of the fluids. The technique of using tracer particles in an Eulerian 
scheme has been used in simulations of core-collapse supernovae [5 6| and Type la supernovae 
0. The initial spatial distribution of the particles follows the mass distribution. The tracer 
particles were evolved by a second-order predictor-corrector method [8|, with the velocity at 
the particle position computed by parabolic interpolation. 

2.1. Verification 

We have partially verified the tracer particle module on a smooth hydrodynamic problem, 
the isentropic vortex. This test problem has a known analytic solution and is suitable for solu- 
tion verification studies (see @, and references therein). We used a uniform grid of 128x128 
equidistant zones and 10,000 randomly distributed tracer particles. In our analysis we focused 
on a subset of three test particles located at r = 0.26, 0.161, and 3.873 cm. We evolved this 
system for t = 100 s, which corresponds to 21 .3 orbital periods for the inner particle and about 
15.4 orbital periods for the middle particle. The outer particle moves very slowly and traveled 
a negligible fraction of its orbit over the simulation. 

Fig.[T](/e// panel) shows the relative error in the particle radius at the final time. The deviation 
from the nominal orbit is generally small and is less than 1% for t < 20 s, or about 4 orbits for 
the inner particle. At later times the error steadily grows but appears to be increasing for only the 
inner particle. The orbit of the inner particle is tight (resolved by about only 20 zones), and the 
general increase in radius likely results from diffusion in the hydrodynamic solver in addition 
to the intrinsic error in tracing the flow. While separating these effects is a future endeavor, 
we are confident that our implementation is correct and that tracer particles closely follow the 
streamlines. 

2.2. A 2-D Deflagration 

The two-dimensional deflagration simulation was performed in cylindrical coordinates and 
azimuthal symmetry on a domain with r = [0,8192] km and z = [-8192,8192] km and an 
effective resolution of 8 km. The initial model was an isothermal white dwarf of temperature 
r wd = 5 x 10 7 K, mass M wd = 1.3 M Q , and radius R w d = 2130 km. The outer boundary condi- 
tions were hydrostatic outflow-only and reflecting boundary conditions were imposed along the 
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Figure 1 . Left panel Relative error in the position of test tracer particles in the isentropic vortex 
verification test problem. Right panel Mass fractions of the isotopes synthesized by the de- 
flagration. Dotted lines connect isotopes of the same element. The number beside each point 
refers to the mass number of that isotope. 



symmetry axis. The deflagration was ignited by making the innermost 100 km the white dwarf 
into completely burned material in hydrostatic equilibrium with its surroundings. The flame 
was then evolved by integrating the hydrodynamic equations supplemented by the advection- 
diffusion-reaction model described above. The deflagration was evolved for 2.0 s, by which 
time « 0.4 M Q had burned. 

3. Nucleosynthesis 

The isotopic abundances were calculated by integrating a reaction network along the trajec- 
tories in (p, T) space followed by the tracers. We use a network consisting of 214 isotopes 
ranging from neutrons, protons, and a-particles to 70 Zn. The coupled ODEs are integrated us- 
ing an adaptive semi-implicit scheme with a sparse matrix solver |[T0"|| . We use the reaction rate 
compilation reaclib 1111 . and incorporate the tabulated weak rates from lfT2l. Screening is com- 
puted according to the prescription of JT3l . Since each particle is independent, the calculation 
is embarrassingly parallel, and large numbers of particles can be quickly post-processed on a 
multi-CPU machine. 

Fig.\J\(right panel) shows the mass fractions of the isotopes synthesized in the deflagration. 
One notes that most of the original C-0 did not burn and that roughly equal masses of 56 Ni 
and 28 Si were made. Only about 0.15 M Q of Ni-peak elements are synthesized. This is much 
less than that found in one-dimensional calculations such as W7 ifbUl and recent 3-D octant 
calculations 0. For our chosen configuration, the burning is incomplete. Large plumes form 
and rise toward the surface of the white dwarf. The large scale circulation arising from these 
flows tends to homogenize the composition as a function of radius: unburned 12 C and 16 0, 
intermediate mass elements and Ni-peak elements are present at roughly constant abundances 
from the center to about 0.9 of the WD mass (as far as the rising plumes had risen for this 
simulation). 
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We did repeat the calculations with a starting abundance that included Ne. From charge and 
mass conservation, one can show that the mass of 56 Ni should decrease linearly with the initial 
mass fraction of 22 Ne [fT51. Calculations with an initial 22 Ne mass fraction of 0.02 and 0.06 also 
find a linear dependence of the mass fraction of 56 Ni, but with a 10% steeper slope. Because the 
underlying hydrodynamical model was the same, however, we cannot draw further conclusions 
about the relationship between ejected 56 Ni mass and progenitor metallicity. We also checked 
how the computed abundances scaled with the number of tracer particles. For a subset of 5000 
randomly selected particles, the difference in the 56 Ni mass from that computed with 10,000 
particles is 8%. 

We finally note that although only about 10% of the star is burnt to Fe-peak elements of that 
fraction, almost all of it is in the form of 56 Ni; the mass fraction of 58 Ni, 54 Fe and 56 Fe are 0.019, 

0. 009. and 0.0005, respectively. This is in contrast to 1-D models, for which electron captures 
convert the innermost « 0.2 M Q to 56 Fe ifTBI . The rising plumes of hot ash carry the burnt 
material to lower density, where electron capture rates are lower. The effect of rising bubbles 
of hot fuel to mix the white dwarf and leave masses of C-0 fuel unburnt has been noted by 
others (for example, [ 17 1) and indeed, is a motivation for suggesting a subsequent transition to 
detonation. 

In conclusion, we have implemented particle tracers in an Eulerian, adaptive mesh hydrody- 
namics code, verified that they closely followed streamlines, and used them for post-processing 
a 2-D simulation of a deflagrating C-0 white dwarf. Future 3-D simulations will use this capa- 
bility for comparison against other simulations [7| and observations. 
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